The General Boson Normal Ordering Problem 
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We solve the boson normal ordering problem for F Uat) r a a "\ , with r, s positive integers, [a, a*] = 
1, i.e. we provide exact and explicit expressions for its normal form TV jf 1 [(a^) r a s ] }, where in 
J\f (F) all a's are to the right. The solution involves integer sequences of numbers which are gener- 
alizations of the conventional Bell and Stirling numbers whose values they assume for r = s = 1. 
A comprehensive theory of such generalized combinatorial numbers is given including closed-form 
expressions (extended Dobinski - type formulas) and generating functions. These last are special 
expectation values in boson coherent states. 



Consider a function F(x) having a Taylor expansion around x = 0, i.e. F{x) — X^feLo F< k\^ x> * • ^ n * n ^ s n0 ^ e we 
will collect the formulas concerning our solution of the normal ordering problem for F [(a^) r a s ] , where a, are the 
boson annihilation and creation operators, [a, = i ; and r and s are positive integers. The normally ordered form 
of the operator F [(a^) r a s ] is denoted by Af {F [(a^a 8 ] }, where in Af(F) all the a's are to the right. It satisfies the 
operator identity: 



Af {F [(a t ) r a s ] } — F [(a t ) r a s "l 



(1) 



Furthermore, an auxiliary symbol : 0(a, a') : will be used, which means expand O in powers of a and and order 
normally assuming they commute [1],[2], 

The combinatorial numbers S(n,k), known as Stirling numbers of the second kind, and their sums B(n), the Bell 
numbers, arise naturally in the normal ordering procedure for r = s = 1, as follows [3]: 



{a)a) n = ^S{n,k){a)) k a k , 



k=l 



B(n) = J2S(n,k), 



(2) 



(3) 



k=l 



which may be taken as definitions of S(n, k) and B(n). In the present work we shall treat the case of general (r, s), 
consequently generalizing these combinatorial numbers. 

For the moment we restrict ourselves to the case r > s, the other alternative being treated later. We do not give 
the proofs of the formulas here; they will be given elsewhere [4]. The case r = s = 1 is known [1],[2] and some features 
of the r > 1, s = 1 case have been published [5]. When needed we shall refer to [1],[2] and [5] where particular cases 
of our more general formulas are discussed. 

We define the set of positive integers S r , s {n, k) entering the expansion: 



s) 



ns 



(4) 
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for n — 1,2,... . Once the S r ^ s (n,k) are known, the normal ordering of [(a^) r a s ]™ is achieved. The same applies to 
any operator F [(a t ) r a s ] if F(x) has a Taylor expansion around x = 0. The row sums of the triangle S r ^ s (n 7 k) are 
given by: 



B r , s (n) = ^2 S r>s (n, fc), 

k—s 

which for any t extend to the polynomials of order ns defined by: 

ns 

B r>a {n,t) =^2S r , s (n,k)t k . 



(5) 



(6) 



All the subsequent formulas are consequences of the following relation, linking the polynomial of Eq.(6) to a certain 
infinite scries: 



oo n 



k=s ' j=l 



...■(k+(j-l)(r-s)-s + l)]t k = 



(7) 



£Sr,.(n,k)t' 



n = 1,2... 



k—s 



which for r, s > 1 and t = 1 is an analogue of the celebrated Dobinski relation, which expresses the combinatorial Bell 
numbers B\.i(n) as a sum of an infinite series [6], [7], [8]: 



, N 1 Af;™ 



(8) 



For a review of various characteristics of Bell numbers B\ t i{n) see [9]. It seems that Bell and Stirling numbers are 
now beginning to appear in textbooks of mathematical physics [10]. 

By setting t — 1 in Eq.(7) we obtain expressions for the generalized Bell numbers B r ^ s (n), n = l, 2... : 



r > s : B, %s (n) = B r , s (n, 1) = 



1 oo n 

Eiill^ + CJ- 1 ^--)! 4 



e f- fc! . 

/c— S J — 1 



— s ^(n-l) 



(r - s) 



E 

fe=0 



fc! 



±± r(1+ fe±i) 

7 — 1 v r — s J 



r = s : B r>r {n) = B r%r (n, 1) = 



e ^ fc! 

fe=0 



(fc + r)! 



fc! 



n-l 



(9) 
(10) 

(11) 



For all r, s we set £? riS (0) = 1 by convention. In Eq.(9) we use the notation m- := m(m — l)(m — 2) . . . (m — s + 1) 
for the falling factorial [11], and in Eq.(10) T(y) is the Euler gamma function. 

The numbers S r ^ s (n, k) are non-zero for s < k < ns, and satisfy by convention S r;S (n, 0) = £„ ; o- We shall refer to 
them as generalized Stirling numbers of the second kind. Their exact expressions in the form of a finite sum are: 



(-1)' 

fc! 

(-1)' 
fc! 



E^n p jUtp+ti- 1 ^- 3 ^ 

p=s ^ j=l 



s-1 



a-^-El „ )(-*)' 

p=0 



(12) 
(13) 



J ' x=l 

from which for r = s = 1 one obtains the standard form of the classical Stirling numbers of the second kind [8] : 

\k k 



Sl!li n,k)= { -jLj2(-iy( k ] P « 

p=l ^ ' 



(14) 
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and for r = 2, s = 1 

Q„ . ^ n — , 

fc! V k-1 

which are the so-called unsigned Lah numbers [5], [8]. 

That the numbers S r , s (n, k) are natural extensions of Si,i(n, k) can be neatly seen by observing their action on the 
space of polynomials generated by falling factorials. Let x be any real number. Then, 

n ns 

]J(x + (j - l)(r - = S rAn, k)xK (16) 

j = l k=s 

This last equation when specified to r = s gives a particularly transparent interpretation of S r ^ r (n, k) as connection 
coefficients: 



(a£) n = 5^S riT .(n,fc)a£. (17) 

A;— r 

For r = s = 1 it boils down to the defining equations for £1,1(71, A;) in terms of x-, used by J. Stirling himself [11]: 

n 

x n = ^S*i,i(n,fc)x^. (18) 
fc=i 

The relations for S r ^ s (n 7 k) become particularly appealing for r — s: Eq.(12) can be reformulated as: 

p—r v 

which differs from Eq.(14) in that the sum in Eq.(19) starts from p = r and p n becomes (p L ) n ■ The recurence relations 
for S rtr (n, k) are the following: 

SV,r(l,»-) = 1, S r>r (n,k) = 0, k<r, nr<k<(n + l)r, (20) 

S r , r (n+l,k)=J2( k+P ~ r \rPS r , r (n,k+ P -r), r<k<nr, n > 1, (21) 

where the definition r— := 1 is made. Note that for r = s = 1 we get the known relation for conventional Stirling 
numbers i.e.: Si,i(n + 1, &) = kSi^(n, k) + Si,i(n, fc — 1) with appropriate initial conditions [11]. 

The "non-diagonal" generalized Bell numbers B r ^ s (n) can always be expressed as special values of generalized 
hypergeometric functions p F q . Algebraic manipulation of Eqs.(9)-(11) yields the following examples. For r > 1, s = 1, 
B r s(n) is a combination of r — 1 different hypergeometric functions of type iF r _i(. . . ;x), each of them evaluated at 
the same value of argument x = (r — l) 1_r ; here are some lowest order cases: 

B2,i («) = ^i^i(r i + l;2;l) = (n-l)!Li 1 2 1 (-l), (22) 

etc. In Eq.(22) Lm\y) is the associated Laguerre polynomial. 

Similarly the series B2 r , r { n ) can be written down in a compact form using the confluent hypergeometric function 
of Kummer: 

B 2r , r (n) = ^- 1 F 1 (rn + l,r + l;l). (24) 
e • r! 

In contrast, the "diagonal" numbers B r ^ r (n) of Eq.(ll), which also can be rewritten as: 

oo 

Br.rH = ~ E TFT _ -i \ | [Kfc + 1) • • • (fe + r - 1)]" , n=l,2,... (25) 
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cannot be expressed through hypergeometric functions. However, the sequences B r ^ r (n) possess a particularity that 
they can always be expressed in terms of conventional Bell numbers and r-nomial (binomial, trinomial. . . ) coefficients. 
For example: 



S 2 , 2 (n)=£( n fc 1 )s 1;1 (n + fc). 

k=0 ^ ' 



(26) 



Some low-order triangles of S r ^ s (n 7 k) and their associated B r ^ s (n) are presented in Table f. 

It turns out that various generating functions for B r ^ s {n) and S r ^ s (n 7 k) can be related to special quantum states, 
which are called coherent states [1], defined as linear combinations of the eigenstates of the number operator, N — at a, 
N\n) — n\n) 7 (n\n') — 5 n ,n' and defined as: 



oo 



\*) 



_i£i± Z " I \ 



(with (z\z) = 1), for z complex. The states \z) satisfy: 

a\z) — z\z). 

It has been noticed in [13] that, 

(z\(Ja) n \z) ^'t 1 Bi,i(n) 

and the exponential generating function (egf) of the numbers Bi^{n), i.e. J2^=o Bis(n)jj satisfies 

n=0 

which is a restatement of the known fact that [1],[2],[5]: 

e Aa+a _ ^Aa 1 ^ _. g a t a(e A -l) . 



00 \ n 

(z\e Xa ' a \z) = (z\ : e «Me A -D : \z) ^'t 1 e^" 1 = £ B hl (n)^ 



How docs it generalize to r, s > 1 ? 
Eqs.(4) and (5) give directly: 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 



valid for all r, s. 

We first treat the case r = s = 1,2.... Note in this context that a Hermitian Hamiltonian (a^) r a r is of great 
importance in quantum optics, as for r = 2, 3, . . . it provides a description of non-linear Kerr-type media[12]. Using 
Eq.(19) we can find the following egf of S r , r (n 7 k): 



n=[k/r] P=r ' ' ' 



(33) 



where \y\ is the ceiling function [11], defined as the nearest integer greater or equal to y. This yields via exchange of 
order of summation: 



e A(at)^ = N ( e A(atr„^ = 1+ : ( ^ ) {e^^"^^ - 1 

fe=r ' .p=r \P / 



(a t a) fc : 



and 



(z\e^\z) |2 L =1 1 + E ^ E(-D P ( J ) K - 



k—r 



(34) 



(35) 
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Another case which can be written in a closed form is r > 1, s = 1, for which the egf for 5 r> i(n, k) reads: 

00 n 1 

E ^ r , 1 (»,Jfe) = ^{(l-(r-l) af )-^-l}. 

n— [fc/r] 

One then obtains [5], [14]: 

e A(at)"a = ^ ( e A(a+)"a^ = . exp j [(J _ A(a t ) r " 1 (r - 1))~^ - l] a t a | : 

For arbitrary r > s the following formula is still valid: 

00 \n / ns \ 

e^ aS = N (e x W) =l+:J2 E S ^ k)(a^a) k : 

n=l \k=s I 

where the explicit form of S r , s (n, k) may be used, see Eqs.(12)-(13). 
A close look at Eqs.(35) and (32) reveals that 

n=0 



(36) 



(37) 



(38) 



(39) 



Comparing Eqs.(39),(37) and (35) we conclude that for r > 1, s = 1 and r — s the egf of respective generalized Bell 
numbers are special matrix elements of e A ( at ) r<lS in coherent states. 

However, for arbitrary r > s > 1 examination of Eqs.(9)-(10) confirms that the numbers B r ^ s (n) increase so rapidly 
with n that one cannot define their egf's meaningfully. In particular, this is reflected by Eq.(38) which is true in its 
operator form, but one needs to exercise care when calculating its matrix elements as the convergence of the result 
may require limitations on A. 

A well-defined and convergent procedure for such sequences is to consider what we call hypergeometric generating 
functions (hgf), i.e. the egf for the ratios S r;S /(n!)*, where t is an appropriately chosen integer. A case in point is 
the series B 3y2 (n) which may be written explicitly from Eq.(lO) as: 



oc 

B 3 ,2W = -E 



(n + k)\(n + k+ 1)! 
Q fc!(fc + !)!(* + 2)! ' 



Its hgf G3.2W is then: 



0ttW _ f I g ^ 2 W + 2, fc + 1; 1; A). 



n=0 



(40) 



(41) 



Similarly one can work out other examples of hgf. See [15] for other instances where this type of hgf appear. 

All the previous considerations also apply to the case r < s. In this case we define the generalized Stirling numbers 
by (note the difference with Eq.(4)): 



r < s : 



[{a*) r a'] n =Af{[(a^) r a s ] n } 



_k—r 



a 



n(s—r) 



By taking the Hcrmitian conjugate of Eq.(4), with the change r <-» s, we obtain 



r < s : 



[(at)^ s ]" 



,k—r 



a 



n(s — r) 



whence the symmetry relation: 

S rtS (n, k) = S Syr (n, k), 



s < k < ns, (r < s). 



(42) 



(43) 



(44) 



The normal ordering of F[a s (a') r ] involves the antinormal to normal order, in the terminology of [16]. With this 
in mind we define the anii-Stirling numbers of the second kind S r ^ s (n, k) as: 



r > s : 



[a s (a r ) r } n = (at)"^- 8 ) 



5^S a , r (n,fc)(at) fc a* 



.k=0 



(45) 
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and similarly aniz-Bell numbers: B r s {n) — Y^kLo S r ,s(n,k) for r > s. 

It can be demonstrated [4] that they satisfy the following symmetry properties: 



S s , r {n, k) — S r , s (n, k), < k < ns, (r > s), 

S r . s (n, k) = S riS (n + 1, k + s), < k < ns, (r > s), 



(46) 
(47) 



and B rtS {n) = B r ^ s (n + 1). 

If in Eq.(l) F(x) has a Taylor expansion of F(x) around x ^ 0, then 




k=0 



(48) 



and our above results should be supplemented by the appropriate binomial expansion coefficients. 

All the B ryS (n) described here are the n-th moments of positive functions on the positive half axis; some solutions 
of the associated Sticltjcs moment problem have been recently discussed at some length [17]. In addition we have 
found that even larger classes of combinatorial sequences are solutions of the moment problem and they can be used 
for the construction of new families of coherent states [18]. 

We are considering the possible combinatorial interpretation of the above results. 
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Table 1. 

Triangles of generalized Stirling and Bell numbers, as defined by Eqs.(4) and(5). 



r = l,s = 1 r = 2,s = 1 
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i(n, k), 1 < k < n 


Bi,i(n) 






52,i (n, k), 1 < k < n 




Sa,i( 
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= 1 
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= 1 
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n 


= 2 


1 


1 


2 


n 


= 2 


2 1 






3 


n 


= 3 


1 


3 1 


5 


n 


= 3 


6 6 


1 




13 


n 


= 4 


1 


7 6 1 


15 


n 


= 4 


24 36 


12 1 




73 


n 


= 5 


1 


15 25 10 1 


52 


n 


= 5 


120 240 


120 20 


1 


501 


n 


= 6 


1 


31 90 65 15 1 


203 


n 


= 6 


720 1800 


1200 300 


30 1 


4051 



r=2, s=2 



S 2 , 2 (n,k), 2 < k < 2n B 2 , 2 (n) 
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= 1 
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= 2 
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1 
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= 3 
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32 


38 


12 
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87 
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= 4 
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208 


652 


576 


188 


24 


1 






1657 


n 


= 5 


16 


1280 


9080 


16944 


12052 


3840 


580 


40 


1 


43833 


n 


= 6 


32 


7744 


116656 


412800 


540080 


322848 


98292 


16000 


1390 60 1 


1515903 



r=3, s=2 

S 3 , 2 {n,k), 2 < k < 2n B 3 , 2 (n) 
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= 1 
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= 2 


6 


6 


1 














13 


n 


= 3 


72 


168 


96 


18 


1 










355 


n 


= 4 


1440 


5760 


6120 


2520 


456 


36 


1 






16333 


n 


= 5 


43200 


259200 


424800 


285120 


92520 


15600 


1380 


60 


1 


1121881 


n 


= 6 


1814400 


15120000 


34776000 


33566400 


16304400 


4379760 


682200 


62400 


3270 90 1 


106708921 



